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Conv ersion  Factors,  U.  S ._  Cu stomary  to  Met r i c  (SI) 
llr Its  of  Measurement 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  converted  to 


metric  (SI)  units  as  follows: 

Multiply 

To  Obtain 

feet 

0.3048 

metres 

inches 

2.54 

centimetres 

kips  (force)  per  square  inch 

6.894757 

megapascals 

knots  (international) 

0.5144444 

metres  per  second 

pounds  (force)  per  square  inch 

6894.757 

pascals 
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PREDICTION 


KjI 


PAVEMENT  ROUGHNESS 


Introduction 

For  many  years  the  emphasis  in  pavement  research  has  been  directed 
toward  the  development  of  a  mechanistic  approach  t.  pavement  design.  Nearly 
all  such  design  procedures  have  been  based  on  an  empirical  correlation  between 
pavement  distress  and  computed  pavement  response  parameters.  The  computed 
response  parameters  are  normally  based  on  a  deterministic  evaluation  of  mate¬ 
rial  properties  and  pavement  section  geometries.  These  design  procedures  do 
represent  an  advancement  in  the  state  of  the  art  of  pavement  design  in  that 
they  provide  a  mechanism  for  considering  basic  material  properties  and  pave¬ 
ment  response.  A  shortcoming  of  these  procedures  is  the  fact  the  procedures 
do  not  lead  directly  to  a  prediction  of  the  functional  performance  of  a 
pavement.  For  military  airfields  the  concept  of  functional  requirements  is 
particularly  important.  The  importunes  of  functional  requirements  are  illu¬ 
strated  in  the  studies  of  methods  for  repairing  aorob  damage  and  the  minimum 
pavement  requirement  for  alternate  pavements  to  be  used  only  in  emergency. 

In  each  of  these  examples  the  conventional  definition  of  pavement  distress 
becomes  meaningless.  The  principal  consideration  is  the  effect  the  pavement 
surface  has  on  the  operating  aircraft.  An  example  of  the  acceleration  forces 
generated  bv  a  B-S2  aircraft  operating  on  a  hypothetical  sinusoidal  profile  is 
shown  in  Fig.  1  (Horn  1977).  For  the  hypothetical  sinusoidal  profile  the 
critical  wave  length  for  various  aircraft  speeds  is  shown  in  Fig.  2  (Horn 
1977) .  Thus  it  is  seen  that  the  forces  generated  are  dependent  on  both  the 
amplitude,  frequency  of  rutting,  and  speed  of  the  aircraft  With  this  in  mind 
then  the  prediction  of  pavement  roughness  is  seen  to  involve  not  only  the 
prediction  of  a  rut  depth  but  also  the  prediction  of  the  distribution  of  the 
ru  1 1 ing . 

Scope 

This  report  presents  a  general  methodology  for  predicting  pavement 
roughness  and  a  specific  methodology  for  predicting  rutting  in  flexible 
pavements.  No  specific  method  was  developed  for  predicting  of  roughness  for 
rigid  pavements  nor  for  the  roughness  caused  by  swell  and  frost  heave  of 
pavements . 


The  prediction  of  pavement  roughness  involves  the  generation  ol  hypo¬ 
thetical  pavements  that  meet  known  statistical  parameters  for  suriace  profile, 
thickness  and  material  properties.  For  each  of  these  pavements  a  mechanistic 
model  is  used  to  predict  the  development  and  distribution  of  rutting  with 
traffic.  Thus  at  any  time  in  traffic  a  profile  for  each  pavement  can  be  gen¬ 
erated  and  the  resulting  profile  analyzed  as  to  the  effect  on  using  aircraft 
by  using  an  aircraft  simulation  program,  such  as  the  one  described  by  Horn 
(1977).  If  a  sufficient  number  of  hypothetical  pavements  are  analyzed,  then  a 
distribution  of  roughness  parameters  can  be  generated  such  that  the  probability 
that  a  certain  roughness  develops  is  predicted. 

Generation  of  Hypothetical  Pavements 

The  generation  of  a  series  of  parametric  values  that  meet  statistical 
data  may  be  n.vompl ished  in  a  number  of  ways.  One  of  the  simplest  is  the  use 
of  a  distribution  function  bv  which  a  series  of  values  will  be  distributed. 

The  values  can  be  selected  randomly  from  the  distribution  such  that  resulting 
set  will  have  no  order  but  will  have  a  specified  mean  and  standard  deviation. 

The  most  commonly  used  function,  the  normal  distribution,  has  some  undesirable 
properties,  mainly  that  the  possible  values  are  unbounded  and  the  data  are 
symmetrically  distributed.  Possibly  a  better  distribution  function  would  be 
the  Beta  function  as  described  by  Harr  (1977).  The  Beta  function  overcomes 
some  of  the  disadvantages  of  the  normal  distribution  in  that  it  may  be  skewed 
and  has  a  specified  minimum  and  maximum  value.  A  computer  program  was  written 
(Appendix  A)  that  provides  a  set  of  values  that  are  distributed  according  to 
the  Beta  function  and  has  a  given  mean,  standard  deviation,  minimum  value  and 
a  maximum  value. 

Consider  lor  example  the  generation  ot  40  values  having  a  mean  of  f> .  7 , 
a  standard  deviation  of  2.2,  a  mini  mum  value  of  2  and  a  maximum  value  of  If. 

The  values  generated  by  the  computer  program  having  the  required  statistical 
properties  are  given  in  Appendix  A.  The  distribution  of  the  values  is  shown 
in  Fig.  1.  Such  a  distribution  could  verv  wall  represent  the  values  of  soil 
strength  in  terms  of  California  Rearing  Ratio  (CUR)  where  CBR  values  are  bounded 
and  have  a  skewed  distribution.  The  parameters  describing  a  pavement  section 
can  also  be  generated  using  the  Beta  distribution  him  lion. 


From  .1  statistical  point  the  randomly  general  od  values  satisfy  t  ho  re¬ 
quirement  of  the  problem.  If  in  the  generation  of  the  properties  for  a  hypo¬ 
thetical  pavement,  the  points  for  which  properties  are  being  generated  are 
sufficient  Iv  far  apart  that  the  properties  are  independent  of  one  another,  then 
the  procedure  is  correct.  For  airfield  pavements  where  the  critical  wavelength 
is  rather  long  and  properties  are  on  1 v  generated  for  these  points, then  the 
routine  mav  be  sufficiently  correct  to  use.  If  the  shorter  wavelengths  are 
to  be  considered  and  the  properties  at  one  point  are  influenced  bv  the  propertie 
of  adjacent  points, then  additional  restraints  must  be  placed  on  the  generation 
scheme . 

Dr.  Per  I'llidtz  developed  such  a  routine  for  prediction  of  roughness  in 
highway  pavements  where  the  critical  wavelengths  are  much  shorter  than  in 
airfield  pavements.  Dr.  i'llidtz  first  tried  what  he  called  a  random  walk  in 
which  the  mean  of  a  parameter  distribution  was  set  equal  to  the  value  ol  the 
previous  point.  This  procedure  resulted  in  unrealistically  smooth  pavements 
and  thus  the  scheme  was  modified  to  what  Dr.  Ullidtz.  named  the  modified  random 
walk.  This  scheme  is  different  in  that  the  mean  of  the  parameter  distribution 
is  a  projection  of  the  values  of  the  previously  determined  two  values.  The 
two  schemes  are  illustrated  in  Tig.  4  (Fig.  2  of  Ullidtz). 

In  a  study  ol  the  significant  characteristics  of  runwav  roughness,  Berens 
and  Newman  (Berens,  1971,  AFFDI.-TR— 73-109)  developed  a  computer  program  tor 
generating  hypothetical  runway  profiles  that  had  a  specified  power  spectral 
density  (PSD). 

also  the  PSD  of  a  number  of  different  pavements  were  determined  and  a 
typical  new  pavement  PSD  was  estimated.  Thus  using  the  computer  program, a 
profile  could  be  determined  that  would  simulate  a  new  pavement.  In  the  simula¬ 
tion  process  a  methodology  was  developed  for  random  spacing  of  the  different 
waves  along  the  pavement  profile.  An  example  of  the  results  of  the  simulation 
process  is  shown  in  Fig.  5. 

The  problem  of  developing  probabilistic  information  about  the  engineering 
parameters  of  a  space  also  exists  in  other  areas  of  geotechnical  engineering. 
Consider  the  paper  by  Wu  and  Wong  (Wu  and  Wong  1981)  that  describes 
the  case  history  of  a  problem  in  prohab i I ist ics  soil  exploration.  In  this 
ease,  the  soil  properties  are  measured  at  several  locations.  The  probability 


contours  ore  developed  based  on  a  proximity  rule.  For  pavement  design,  it  is 
usually  not  possible  to  have  measured  parameters  since  the  pavement  system  is  to 
be  constructed  in  the  future.  Still  such  a  proximity  rule  could  be  used  partic¬ 
ularly  if  short  wavelengths  are  important.  The  basic  concept  would  be  to 
randomly  select  parameters  at  distances  such  that  they  can  be  considered 
mutually  independent.  These  values  can  then  be  considered  the  same  as  measured 
values  and  a  proximity  rule  used  to  develop  parameters  between  the  randomly 
determined  values. 

Prediction  of  Pavement  Rutting 

Genera  1 

Two  basic  approaches  were  available  for  prediction  of  the  pavement 
rutting  -  the  statistical  and  mechanistic.  Barber  et  al.  (1978)  conducted  a 
study  of  numerous  pavements  and  developed  a  statistical  model  for  prediction 
of  the  rut  depth  as  a  function  of  CBR,  thickness,  and  traffic.  This  particular 
model  provides  the  rut  depth  as  a  deterministic  value  or  a  statistical  value. 

If  the  rut  depth  is  determined  statistically,  i.e.  in  terms  of  a  mean  and  variance, 
then  the  rut  depth  can  be  distributed  directly.  Although  Barber's  model  cur¬ 
rently  exists,  for  reasons  of  verr.itility  and  laboratory  testing  considerations 
a  mechanistic  model  for  prediction  of  rutting  was  developed. 

Permanent  Strain  Mode  1 

The  mechanistic  model  is  based  on  a  laboratory-determined  relationship 

between  the  permanent  deformation  (e  )  and  the  resilient  deformation  (t  )  and 

P  k 

state-of-stress  along  a  vertical  axis  of  the  pavement  system.  The  first  work 
accomplished  toward  developing  the  mechanistic  model  was  to  review  available 
laboratory  data  for  the  development  of  the  relationship  between  e  and  and 
to  show  the  feasibility  of  using  the  approach  for  predicting  rut  depth.  This 
study  was  rather  lengthy  and  involved;  therefore  the  results  are  presented  sepa¬ 
rately  In  Appendix  B.  Although  the  methodology  presented  in  Appendix  B  was  a 
viable  methodology,  it  did  not  lend  itself  to  a  computerized  procedure. 

From  the  data  analyzed  in  Appendix  B  it  was  apparent  that  the  relationship 

between  e  and  £  is  influenced  by  the  soil  strength.  To  develop  the  relation- 
P  R 

ship  with  soil  strength  the  WES  data  were  plotted  as  shown  in  Fig.  6.  For  each 

of  the  repetition  levels  a  relationship  is  shown  between  the  ratio  of  e  and  c 

P  K 

and  the  resilient  modulus  (1L) .  By  considering  the  relationship  shown  in  Fig.  6, 
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che  relationship  that 
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where 

R  =  0.4  (Stress  Reps)^'^" 


a  =  repeated  deviator  stress  in  laboratory 
triaxial  test 

t  =  measured  resilient  strain  in  laboratory 
triaxial  test 


is  developed. 

As  an  independent  check  of  the  relationship,  a  comparison  was  made  with 
experimental  data  reported  by  Ogawa  (19721.  This  comparison  is  shown  in  Fig.  7. 
Model  for  Resilient  Modulus 

For  the  purpose  of  this  study  a  comparison  with  actual  test  data  would 
be  most  desirable.  One  problem  is  that  the  strain  model  is  based  on  where 
the  measured  material  property  in  the  field  test  sections  is  the  CBR.  Making 
the  conversion  from  CBR  to  Mr  is  rather  difficult  since  the  is  a  function 
of  stress  as  well  as  material  properties.  This  conversion  was  made  by  devel¬ 
oping  an  empirical  correlation  between  field  CBR  and  M^.  The  correlation  may 
be  expressed  bv  the  equation 

\  =  10* 

where 

*  ■  4-5682  -  sHrlk  «"5  *  V 

a,  =  repeated  deviator  stress 
d 


The  relationship  is  shown  graphically  in  Fig.  8. 

This  correlation  was  developed  from  data  obtained  in  connection  with  a 
study  (Farker  et  al.  1979)  to  develop  a  correlation  between  plate  bearing  value 
and  M^.  For  the  study  a  number  of  field  sites  were  selected  based  on  an  attempt 
at  having  sites  covering  a  wide  range  of  soil  types.  At  each  site  plate  bearing 
and  CBR  tests  were  conducted  and  undisturbed  soil  samples  were  taken.  Resilient 
modulus  tests  were  conducted  by  the  WES  Geotechnical  Laboratory.  An  example  of  the 


Idbwr.uurv  tt-.si  1  lent  modulus  data  and  roni|>.ir  ison  with  the  ninth- 1  is  given  in 
Fi«.  9.  U'itli  t lie  models  now  it  is  possible,  provided  the  stress  distribution 
is  known,  to  make  prediction.-;  of  permanent  deformation  for  various  prototype 
pavement  tests  that  have  lit  on  conducted  at  WFS . 

Stj-ess  Mode  1 

Mien  considering  the  stress  distribution  in  flexible  pavements,  the 
structural  layers  above  the  subgrade  are  divided  into  bound  and  unbound  materials. 
Bound  materials  are  capable  of  substaining  tensile  stress  and  thus  tend  to  dis¬ 
tribute  load  as  a  slab;  however,  complete  slab  action  is  unlikely.  Shrinkage 
and  load  associated  cracking  reduce  the  load-distributing  characteristics  of 
stabilized  layers  in  pavement  systems.  In  another  study  (the  report  to  be 
published  later)  it  was  shown  from  data  obtained  in  prototype  testing  of  air¬ 
field  pavements  that  bound  layers  give  a  load-distributing  effect  equivalent  to 
a  two-layered  system  having  a  modulus  ratio  (modulus  of  elasticity  of  the  top 
layer  divided  by  the  modulus  of  elasticity  of  the  bottom  layer)  between  3  and 
4.  This  ratio  will  depend  somewhat  on  the  type  material  being  stabilized  and 
the  amount  of  stabilization.  For  graveLs  that  are  well  stabilized  the  ratio 
will  approach  4;  whereas  for  clays  the  ratio  will  be  closer  to  3.  If  there  is 
very  dense  cracking  in  the  stabilized  material,  the  ratio  could  be  lower  but 
this  is  not  likely  to  occur  before  complete  failure  of  the  pavement  system. 

Unbound  structural  layers,  both  crushed  and  uncrushed  gravels,  in  a  pave¬ 
ment  system  distributed  the  load  equivalent  to  a  two-lavered  system  having  a 
modular  ratio  of  1;  that  is  to  say  the  stresses  will  be  distributed  according 
to  the  Boussinesq  stress  model. 

Since  the  procedure  for  predicting  roughness  requires  the  computation  of 
stress  many  times,  it  is  desirable  to  have  a  rapid  means  for  performing  these 
computations.  For  this  a  computer  program  was  written  using  tabulated  stress 
factors  that  are  a  function  of  the  r/t  ratio,  r/z  ratio,  and  the  modulus  ratio 
where  r  is  the  radius  o’  the  load  area,  t  is  the  thickness  of  the  structural 
layer,  and  z  is  depth  to  the  point  for  computation  of  stress.  An  example  of 
the  stress  distribution  for  the  G-SA  is  given  in  Fig.  10. 

Rutting  Model 

The  rutting  model  consists  of  the  permanent  strain  model  and  the  stress 
distribution  model.  This  model  Is  illustrated  in  Fig.  11.  There  are  some 
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!.•  sumpt  ions  tii.it  must  he  included  with  the  use  of  the  model.  She  first  assump¬ 
tion  is  that  a  stress  repetition  is  synonymous  with  a  coverage.  The  second 
assumption  is  that  no  permanent  strain  occurs  in  hound  layers.  lor  cement  or 
1 ime-s i ah i I i zed  material,  this  has  been  found  to  essentially  be  true  as  lone  as 
the  hound  materials  completely  fail.  ror  asphalt  hound  materials  rutting  can 
occur  within  these  Livers.  this  is  particularly  true  if  a  poor  design  is  used 
in  the  mix.  This  report  presents  no  rutting  model  for  use  in  predicting  rutting 
in  the  asphalt  bound  materials.  Fortunately  most  military  airfields  have  well- 
designed  asphalt  hound  layers  that  are  fairly  thin  and  thus  the  rutting  within 
these  layers  will  be  negligible.  The  third  assumption  is  that  although  the 
permanent  strain  model  was  developed  based  on  testing  of  subgrade  soils,  it  is 
assumed  to  be  useable  for  granular  materials.  This  assumption  can  be  justified 
from  the  data  presented  in  Appendix  B. 

The  rutting  model  as  presented  in  Fig.  11  provides  a  deterministic 
evaluation  of  the  rut  depth.  If  the  model  described  for  generation  of  a  hypo¬ 
thetical  pavement  is  used  for  generation  of  the  material  properties  and  struc¬ 
tural  layer  thicknesses  and  a  sufficient  number  of  sections  are  generated,  then 
the  variations  in  rut  depth  can  he  simulated.  This  procedure,  referred  to  as  the 
Monte  Carlo  procedure,  can  be  used  to  determine  the  probability  of  the  rut  exceed 
ing  a  certain  value  or,  if  used  in  connection  with  an  aircraft  simulation  program 
the  probability  of  the  acceleration  forces  exceeding  a  certain  value. 

Prof.  M.  Harr  introduced  in  a  short  course  a  procedure  he  referred  to 
as  Rosonb 1  net h ' s  procedure  chat  could  he  used  to  compute  the  rut  depth  in  sta¬ 
tistical  terms.  I  ho  procedures  uses  a  finite  difference  procedure  where  the 
dependent  value  of  a  function  is  . ompufed  at  both  the  moan  plus  a  standard 
deviation  and  tire  mean  minus  a  standard  deviation  of  the  independent  variable. 
These  two  values  arc  used  to  compute  the  mean  and  deviation  of  the  d< pendent 
variable.  The  procedure  for  a  single  independent  variable  is  for 

V  =  f (X) 

then  Y  =  1/2  (X  +  X_) 

and  S  -  1/2  (X+  +  X  ) 
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where 


Y  =  mean  of  Y 

S  =  standard  deviation  of  Y 
v 

Y  »  F(X  +  S  ) 

+  x 

Y  =  F (X  +  S  ) 

x 

X  =  mean  of  X 

S  =  standard  deviation  of  X 
x 

The  variance,  Y  ,  of  Y  can  be  computed  by 

S 

V  =  1  (100) 
>'  Y 


This  procedure  can  be  extended  to  multiple  independent  variables  by  considering 
the  functional  value  for  values  of  the  mean  plus  or  minus  a  standard  deviation 
for  each  of  the  independent  variables,  i.e.. 


y 


1 


S 


1 


2  (X!  +  X~) 

\  (x!  +  xl) 


where 

=  the  mean  plus  a  deviation  for  variable  1 
=  the  mean  minus  a  deviation  for  variable  1 

2 

X+  =  the  mean  plus  a  deviation  for  variable  2 
2 

X  =  the  mean  minus  a  deviation  for  variable  2 

jj  =  a  mean  of  the  function  with  respect  to  variable  1  with  all  other 
variables  at  the  mean 

=  a  standard  deviation  with  respect  to  variable  2  with  all  other 
variables  at  the  mean 

v2  =  same  as  y^  but  with  respect  to  variable  2 
S0  =  same  as  but  with  respect  to  variable  2 


Also  the  value  of  the  function  at  the  mean  of  all  variables  is  used  and  is 
computed  by 
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.1*  =  I 

wlie  re 

X^  =  the  mean  of  variable  1 

X~  =  the  mean  of  variable  - 
Now  the  equation 

u  l'  | 
l,*  U*  u* 

i  an  be  used  to  compute  the  mean,  i. ,  ol  the  tun, lion.  For  v.  irian,  ••  t  he  equations 

V  - 

’  i 

V  .  K 

0 

/■ 

are  used  to  compute  the  variance,  Vj  and  V,,,  witli  respect  to  each  of  the  variables 
Now  the  equation 

1  -  V“  =  (1  +  V  ~)  (l  +  V  -’) 

1 

can  be  used  to  compute  the  variance  ot  the  function  with  respect  to  all  variables. 
The  Rosenblueth  procedure  was  combined  with  the  ruttinc  model  and  computerized 
to  give  a  computer  program  for  computing  the  rut  depth  in  statistical  terms. 

The  listing  of  this  program  is  given  in  Apnemlix  C.  With  the  rut  depth  compu¬ 
ted  in  statistical  terms,  the  rut  could  be  distributed  directly  to  produce  the 
runway  profile. 

Other  Cons  ide^jr ci_t:  i_ons 

A  methodology  has  been  presented  for  generating  a  new  airfield  pave¬ 
ment  profile  and  material  properties  that  meet  given  statistical  data.  For 
flexible  airfield  pavements  these  data  can  be  used  to  predict  the  development 
oi  permanent  deformation  with  traffic  and  thus  construct  a  prediction  ot  the 
profile  with  traffic.  Using  the  profile  as  input  in  an  aircraft  simulation 
program,  a  measure  of  roughness  is  obtained.  The  methodology  developed  thus  far 
is  rather  crude  but  does  illustrate  the  feasibility  of  the  approach  and  provides 
basic  models  necessary  for  future  dove! opment .  As  mentioned  earlier,  a  pressing 
need  is  for  the  development  of  the  proximity  rules  to  bo  used  in  generating  the 
initial  pavement  sections.  Verification  of  the  proximity  rules  should  be 
accomplished  by  comparison  with  actual  measured  airfield  data. 
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Also  a  factor  that  could  be  considered  is  the  development  of  struc¬ 
tural  cracking.  Models  are  currently  available  for  both  rigid  and  flexible 
pavements  for  the  prediction  of  structural  cracking,  but  the  effect  of  the 
cracking  on  the  roughness  of  a  pavement  is  a  relatively  unknown  quantity. 

Due  to  the  interaction  of  the  aircraft  and  runway  surface  the  loading 
to  the  pavement  is  not  necessarily  equal  ro  the  stated  loading  of  the  aircraft. 
The  consideration  of  dynamic  load  in  the  procedure  would  add  an  order  of  magni¬ 
tude  in  the  solution  process.  Such  loadings  are  a  function  of  both  profile 
and  speed, and  thus  consideration  of  the  determination  of  the  dynamic  loading 
would  require  continuous  updating  as  the  profile  changed. 

The  profile  could  also  be  affected  by  climatic  conditions,  particularly 
in  frost  areas.  For  pavements  that  are  not  designed  for  full  frost  protection 
the  frost  heave  is  likely  to  be  very  critical.  The  weakening  of  the  subgrade 
during  spring  thaw  is  another  factor  that  should  be  considered  as  the  permanent 
deformation  that  occurs  during  this  period  of  time  may  be  greater  than  the 
deformation  occurring  during  the  remainder  of  the  year. 

For  rigid  pavements  the  joints  pose  a  special  problem.  There  is  some 
statistical  chance  of  dowel  failures,  spalling,  pumping,  or  other  joint  problems 
that  could  create  a  rough  pavement.  For  rigid  pavements,  more  so  than  flexible 
pavements,  research  aimed  at  solving  the  problems  of  predicting  pavement  rough¬ 
ness  is  almost  nonexistent. 

Summary 

A  simplified  procedure  has  been  developed  for  considering  roughness  in 
pavement  design.  The  procedure  used  the  variation  in  pavement  section  Co  gen¬ 
erate  a  stochastic  airfield  pavement.  A  model  was  developed  for  predicting 
the  permanent  surface  deformation  with  aircraft  traffic.  Thus  the  predicted 
deformation  can  be  added  to  the  initial  profile  to  yield  a  pavement  surface  at 
any  point  in  time.  The  predicted  pavement  section  can  be  tested  for  roughness 
by  the  use  of  an  aircraft  simulation  computer  program. 

Extension  of  this  simplified  procedure  to  provide  consideration  of  some 
of  the  other  factors  that  affect  pavement  roughness  would  require  a  large 
research  effort.  Although  the  research  effort  for  complete  development  of  the 
methodology  is  large,  it  is  a  project  that  should  be  undertaken.  The  return 
on  such  an  effort  in  terms  of  better  pavement  design  and  pavement  management 
data  would  justify  the  monies  spent. 
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AIRCRAFT  ACCELERATION 
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Figure  1.  Acceleration  forces  for  hypothetical  sinusoidal 
profile;  B-52H  (after  Horn  1977) 
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critical  pavement  wavelength 


umber  of  Occurrences 


Value 

Figure  3.  Histogram  of  values  for 
Beta  distribution 


Figure  4.  Illustration  of  random  walk  and  modified 
random  walk  (after  Ullidtz  1978) 
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Figure  5.  Illustration  of  generation  of  runway  profile  using 
power  spectral  density  (after  Berens  and  Newman  1973) 
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Figure  6.  Relationship  between  Ep/E^  an d  (Sheet  1  of  3) 
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Figure  6.  (Sheet  2  of  3) 
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Figure  6.  (Sheet  3  of  3) 
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Comparison  of  permanent  strain  model  with  test  data 
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Relationship  for  predicting 


RESILIENT  MODULUS,  MR  (KSl) 


Figure  9.  Example  resilient  modulus  data  and 
comparison  with  model 
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Figure  10.  Stress  distribution  for  C-5A  for 
15-in.  structural  layer 
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Appendix  B:  Development  of  a 
Basic  Rutting  Model 


TNTRGD'JCTier: 

The  Waterways  Experiment  Station  (WE3),  in  seeking  a  more  rational  approach 

to  design  of  flexible  pavements  than  the  presently  used  CBR  design  procedure, 

has  developed,  under  sponsorship  of  the  Federal  Aviation  Administration  (FAA) 

ar.d  the  Office,  thief  of  Engineers  (OCE),  subgrade  strain  criteria  for  both 

roads  ^Brabston,  et  al.  )  and  airports  (Barker  and  Brabston  ).  The  criteria 

now  beiv  incorporated  into  both  FAA  and  CCE  design  procedures  are  considered 

to  be  the  basis  for  limiting  the  rutting  of  flexible  pavement.  The  approach 

of  limiting  the  resilient  strain  in  the  subgrade  in  order  to  limit  rutting 

of  the  subgrade  implies  that  a  relationship  exists  between  the  permanent 

strain  and  the  resilient  strain  of  subgrade  soils.  In  both  of  the  design 

procedures  the  limiting  subgrade  strain  criteria  were  presented  as  a  function 

o 

of  the  subgrade  modulus.  The  author*"  discussed  this  point  ir.  detail  and 
presented  the  results  ox'  laboratory  tests  which  indicated  that  the  relation¬ 
ship  between  the  permanent  strain  and  resilient  strain  was  indeed  a  function 
of  the  stiffness  of  the  material.  The  work  in  developing  these  strain  cri¬ 
teria  and  in  studying  the  results  of  laboratory  tests  for  conformation  of  the 
criteria  created  ati  interest  in  the  relationship  between  the  permanent  strain 
and  resilient  strain  of  subgrade  soils. 

SUBGRADE  STRAIN  CRITERIA 

A  comparison  of  the  different  subgrade  criteria  is  presented  in  Figure  1. 
Ail  of  the  criteria  presented  with  the  exception  of  that  developed  by  Dr.  Chou 
! Reference  d!  has  been  developed  from  pavement  sections  conforming  to  some 
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provi  yi-ly  established  design  standards.  Although  different  design  standards 
r.ay  have  been  used  and  different  methods  were  employed  in  characterizing  the 
granular  materials,  it  is  seen  when  grouped  that  all  of  the  criteria  form  a 
relatively  narrow  band  across  a  wide  range  of  repetition  levels.  Consider 
mat  for  repetition  level  of  1  x  10'’  repetitions  the  range  of  the  criteria  is 
from  0.'.  t;  x  10  "  to  1.!*  x  10  in. /in.  Even  if  the  criteria  of  Finn,  et  aL, 
ar.d  C?u.u  were  extrapolated  to  this  level  of  repetition,  the  range  would  not 
be  increased.  From  the  evidence  it  would  appear  that  for  a  given  level  of 
repetitions  the  resilient  strain  at  which  the  permanent  strain  becomes  unaccept¬ 
able  wou.lu  be  within  a  fairly  narrow  band.  The  data  presented  in  Figure  1 
strongly  indicate  a  unique  relations!’. ip  between  permanent  strain  and  resilient 
strain.  Considering  the  emphasis  being  placed  on  the  use  of  the  repeated  load 
triaxial  test  in  which  both  the  resilient  and  permanent  strains  are  measured, 
it  would  seem  that  the  concept  of  limiting  subgrade  strain  criteria  could  be 
substantiated  or  disputed  from  the  results  of  such  laboratory  tests. 

LABORATORY  TEST 

In  recent  years  much  interest  has  been  generated  principally  by  the  work 

q 

of  Barksdale'  in  the  use  of  the  repeated  load  triaxial  test  as  a  method  to 
predict  rutting  of  a  pavement.  The  aim  of  nearly  all  of  the  laboratory 
rutting  tests  performed  to  date  has  been:  first  ,  to  define  the  permanent 
deformation  as  a  function  of  the  applied  stress;  and  second,  to  define  the 
resilient  modulus  of  the  material  as  a  function  of  stress.  Only  one  experiment 
known  to  the  author  has  been  performed  with  the  first  objective  of  defining 
the  relationship  between  resilient  strain  and  permanent  strain.  This  experiment 
reported  by  Chisolm^  was  conducted  to  substantiate  the  concept  by  the  author 
that  t tie  allowable  resilient  subgrade  strain  is  a  function  of  modulus  of  the 


lient  strain  and 


subgr’>le;  i.e.,  to  establish  the  relationship  between  resi 
permanent  strain  l'or  soils  having  different  moduli.  To  accomplish  this  objec¬ 
tive,  a  Vicksburg  heavy  (buckshot)  clay  (CH  an:  E-ll)  was  molded  at  four  dif¬ 
ferent  water  contents  and  tested  in  a  repeated  load  triaxial  test.  Since  the 
object  of  the  experiment  was  to  determine  only  the  relative  relationships 
between  resilient  strain  and  permanent  strain  between  the  samples,  only  1000 
load  repetitions  were  applied  to  each  sample.  The  results  of  the  experiment 
are  shown  in  Figure  2.  Another  experiment  in  which  the  establishment  of  the 
relationships  between  resilient  and  permanent  strain  was  only  a  minor  part 
was  conducted  at  the  University  of  California  at  Berkeley  by  Ogawa  (Refer¬ 
ence  11).  The  results  of  this  experiment  in  regard  to  the  relationship  between 
resilient  strain  and  permanent  strain  are  shown  in  Figure  3.  The  results  of 
these  tests  are  close  to  those  of  WES,  even  though  the  soils  are  different  and 
the  tests  were  conducted  at  different  strain  levels  and  a  different  number  of 
strain  repetitions  were  applied.  Both  Chisolm  and  Ogawa  used  almost  identical 
equipment  and  procedures  in  conducting  the  two  experiments.  One  seemingly 
important  factor  is  that  both  used  LVDT  holding  clamps  placed  on  the  specimen 
(a  description  of  the  measuring  device  is  given  in  References  10  and  ll)  for 
measuring  both  resilient  and  permanent  strains.  When  examining  the  results 
of  other  researchers, the  evidence  indicates  a  distinct,  difference  in  the  test 
results  depending  on  how  the  strains  are  measured. 

The  majority  of  research  reports  studied  were  written  primarily  to  present 
the  permanent  or  the  resilient  deformation  characteristics  of  a  material  but 
not  the  relationship  between  the  two.  Thus,  in  order  to  develop  these  data, 
in  many  cases  it  was  necessary  to  calculate  missing  parameters  which,  in  most 
cases,  was  the  resilient  strain.  The  data  studied  were  that  reported  by 
Barksdale, ^  Kalcheff  and  Hicks,  ^  Fosnberg,1^  feed  and  McNeill,'1'*  Brown  et  al .  , 
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RESILIENT  STRAIN  eR,  I0‘3IN./IN. 

Figure  B2.  Relationships  between  permanent  strain  and 
resilient  strain  for  buckshot  clay 
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RESILIENT  STRAIN  (r  f  I0'3  IN./IN. 

Relationship  between  permanent  strain  and 
resilient  strain  for  different  materials 


that  the  tost  pcir.tr  i'rr  hi >*hvr  strains  were  ;i  tair.ei  by  insrep.iir.^  the 
water  content  and  vrul  i  therefore  represent  IV. ft  free.  sasp-es  of  lever  streng 
If  the  Soilr  in  thir.  experiment  he  nave  similarly  t :  three  tested  by  Ihisrlm 
.and  0<*ava ,  the  permanent  strain  fer  a  particular  resilient  strain  would  he 
hi. "her  for  the  soil  at  the  hi.ther  water  content  but  lower  stress  than  for  the 
soil  at  the  lower  wafer  content  but  higher  stress.  The  effect  of  the  vary: no 
c.f  the  wafer  content  to  obtain  higher  resilient  strains  would  he  t.  re-  nice 
a  steeper  relationship  than  would  have  been  ob'ained  by  maintaining  a 
constant  water  content  ant  increasing  the  applied  stress  to  obtain  higher 
resilient  strains. 

Another  interesting  but  difficult  set  of  data  to  interpret  was  the 
data  presented  in  Reference  If.  In  this  experiment  ,  the  effect  of  stress 
history  of  the  samples  on  permanent  de format  ion  uni  resilient  strain  were 

B7 


1 

22 

Lucan 

29 

percent. 

22 

percent 

7  1 

2 

Wexford 

percent 

16 

percen4' 

i-  ; 

3 

Coal  Island 

6  7 

percent 

25 

percent 

12  i 

« 

C;l 

Dublin 

35 

percent 

19 

percent 

16 

5 

i'evor.  (TA' 

c  ? 

percent 

23 

percent 

2C 

6 

r*i: 

Darniale  (A'1 

r 

> 

p, or  sent 

27 

percent 

25 

7 

Cl 

Darndale  (E) 

U  6 

percent 

27 

percent 

IS 

8 

CL 

Darniaie  (?) 

39 

percent 

22 

percent 

17 

9 

■J*  L 

Larndale  (0) 

33 

percent 

16 

percent 

2  5 

Soil  type  refers  to  the  Unified  Soil  Classification  System. 


PERMANENT  STRAIN  Ep>  I0'3  IN./IN. 
(AT  10,000  STRAIN  REPETITIONS) 
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RESILIENT  STRAIN  €R  ,  (O3  IN./IN. 


Figure  B4.  Relationships  between  permanent  p+-'in  and 
resilient  from  data  developed  by  nnvran  ot 


B9 


etui  led.  A; thou*:!,  it  was  evrul-.U.-.i  in  the  paper  t:.ut  stress  hist-.ry  had  -.n 
effect  on  the  permanent  strain  t  it  nj  effect  resilient  strain,  the  data 
as  ri-lte:  in  figure  ?  iniicntri  the  relationship  betw-.-en  permanent  strain 
an i  resilient  strain  is  little  affected  tv  n>...  lento  iifferences  in  stress 
history.  In  the  plot,  a  distinct  relationship  for  each  overconsolidat ion 
ratio  is  presented  but  the  difference  appears  i neons: stent  and  could  be  experi¬ 
mental  error.  A  single  relationship  would  probably  suffice  for  this  set  of 
data.  The  data  presented  by  Barksdale  (Figure  3),  feed  (Figure  5),  and 
Fossberg  (Figure  3)  are  all  straight  forward ,  each  of  which  provided  useful 
information.  Of  all  the  tests,  those  conducted  by  Fossberg  covered  the 
largest  range  (.0002  to  .0042  in. /in.)  of  resilient  strains.  Two  tests,  one 
by  Barksdale  and  cne  by  Ogawu,  indicated  very  abrupt  change  in  the  slope  of 
the  relationship.  For  the  test  by  Barksdale  the  abrupt  change  came  at  a 
resilient  strain  of  .0013  in. /in.  and  for  Ogawa  at  .0011  in. /in.  A  mere 
gradual  change  at  about  the  same  magnitude  of  resilient  strain  was  indicated 
in  the  results  presented  by  Fossberg.  Such  behavior  lends  strong  support  to 
the  use  of  the  limiting  subgrade  strain  concept  in  pavement  design.  It  can 
also  be  pointed  out  that  the  strain  criteria  developed  to  date  (presented  in 
Figure  1)  are  consistent  with  the  laboratory  results. 

In  examining  the  data  plots,  it  was  noted  that  results  of  Chisolm, 
Fossberg,  and  Ogawa  group  together  and  the  data  of  the  other  researchers  fall 
into  another  group.  It  has  already  beem  emphasized  that  Chisolm  and  Ogawa 
employed  inside  LVDT  holding  clamps  attached  to  the  specimen  to  measure  both 
resilient  and  permanent  strains.  Fossberg  also  used  the  same  devices  for 
the  tests  he  conducted.  Barksdale  found  that  such  measuring  devices  gave 
him  inconsistent  readings  of  permanent  deformation  and  therefore  he  employed 
outside  measuring  devices  for  measuring  permanent  strain  but  still  used  LVDT 


_r  -  A#.V 
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PERMANENT  STRAIN  ip,  10’ 3  IN./lN 
(AT  1000  STRAIN  REPETITIONS) 


Figure  B5 .  Relationship  between  permanent  strain  and 
resilient  from  data  presented  in  Reference 


Clunra  for  measuring  re:'  i  L  i  <-n‘  strains.  it.  ir.  believed  that  the  other 
researchers  all  use  :  rut:  ’  :<»  me.usur :  rig  d:»vio«>s  which  measured  the  deflection 
over  the  entire  sample.  Thus,  the  data  full  into  two  groups;  that  in  which 
the  permanent  strain  war.  measured  using  LYDT  clamps  and  thus  only  over  the 
center  portion  of  the  s-tnple,  and  that  ir.  which  the  permanent  strain  was 
measured  using  outside  devices  and  thus  measured  the  strain  over  the  entire 
sample.  The  difference  in  the  two  groups  cf  data  is  aline  st  an  order  of 
magnitude . 

Additional  data  are  available  by  including  tests  on  granular  materials. 
Such  data  have  teen  presented  by  Chisolm,  Barksdale,  and  Kalcheff  and  Hicks. 
The  plot  of  there  data  is  shown  in  Figure  6.  These  data  indicate  a  steeper 
relationship  between  resilient  strain  and  permanent  strain  than  was  indicated 
for  subgrade  soils.  Considering  the  great  difference  between  the  material 
properties,  the  relationships  for  the  granular  material  are  surprisingly  close 
to  those  for  subgrade  soils.  If  is  quite  possible  that  the  limiting  strain 
concept  could  be  extended  to  apply  to  granular  subbase  and  base  materials. 

USE  OF  LABORATORY  DATA 

It  has  been  shown  that  the  repeated  lead  triaxial  test  is  a  method  for 
developing  the  relationship  between  resilient  strain  and  permanent  strain. 
There  is  still  some  question  as  to  the  best  procedure  for  measuring  the 
permanent  strain  and  that  the  results  obtained  will  depend  on  the  particular 
procedure  used. 

In  addition  tc  establishing  or  verifying  limiting  strain  criteria  the 
data  presented  in  this  manner  can  a] so  be  used  to  estimate  the  permanent 
deformation  of  the  ;;  upgrade .  In  a  report  on  the  structural  analysis 
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Figure  B6.  Relationships  between  permanent  strain  and 
resilient  strain  for  granular  mate- ials 
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(Reference  l8)  of  insulated  layer:;  the  author  used  the  VKf  lest  data  to  show 


that  the  rutting  of  the  subgradt  was  an  insignificant  part  of  the  total  rutting 
of  tiie  paveneru.  In  the  analysis  the  distribution  of  the  computed  vertical 
resilient  strain  within  the  subgrade  is  shown  in  Figure  7.  Using  this  distribu¬ 
tion  of  resilient  strain  and  the  relationship  between  resilient  strain  and 
permanent  strain,  as  established  by  Chisolm,  then  the  distribution  of  permanent 


strain  as  shown  in  Figure  7  was  determine!  (that  is,  e 


-4  ), 

tR 


where 


e  =  the  computed  permanent  strain 
P 

en  =  the  computed  resilient  strain 
n 


e'  =  measured  permanent  strain  in  repeated  load  triaxial  test 
P 

e'  =  measured  resilient  strain  in  repeated  load  triaxial  test 
The  total  permanent  deformation  (A^)  was  then  determined  by  assuming  the 

permanent  strain  went  to  zero  at  a  depth  of  120  in.  below  the  top  of  the  sub- 

,0 


grade  and  computing  the  area  under  the  cur 


;-ve;  i.e.,  ^  =  j 


R 


120 


E 

_J= 

lcR 


Using 


tiiis  procedure,  the  deformation  at  the  top  of  the  subgrade  was  estimated  to  be 
0.08  in. 

A  conservative  estimate  of  the  permanent  deformation  at  the  top  of  the 

subgrade  can  be  estimated  quickly  rind  easily  by  using  the  ratio  of  permanent 

strain  to  resilient  strain  as  determined  from  the  computed  resilient  strain 

at  the  top  of  the  subgrade.  The  ratio  is  then  used  as  a  constant  multiplier 

to  the  computed  subgrade  resilient  deformation  (A  )  to  estimate  the  subgrade 

n 

permanent  deformation.  This  assumes  that  the  ratio  of  permanent  strain  to 
resilient  strain  remains  constant  with  depth  which,  of  course,  is  not  true. 

E  ' 

r  is  a  constant,  then  the  previous  equation  becomes 
eR 


BlU 


This  is  to  say  if 


E  '  0  £  ' 

A  =  —7  ,  I  £„  =  —7  A„  .  For  the  ex;u:.ile  given  above  the  strain  ratio  is  about 

P  £r  L  r  cr  R 

I. 6  and  the  computed  resilient  deformation  at  the  top  of  the  subgrade  was  0.115, 
giving  an  estimated  permanent  deformation  of  0.l8h.  This  estimate  is  over  twice 
the  previous  estimate.  The  difference  being,  as  mentioned  before,  that  the 
maximum  strain  was  assumed  constant  with  depth  and  also  that  previously  it  was 
assumed  that  the  permanent  deformation  was  zero  at  120  in.,  whereas  in  the 
latter  estimate  it  was  assumed  the  permanent  strain  extends  to  an  infinite 
depth . 

The  same  procedure  can  be  used  to  compute  the  permanent  deformation  in 
other  layers  of  the  pavement  systems.  Consider  the  three  pavement  types  as 
shown  ir.  Figure  8  in  which  sections  1  and  2  were  subjected  to  simulated  traf¬ 
fic  of  a  C-5A  aircraft  and  section  3  to  the  simulated  traffic  of  a  7^*7  aircraft. 

To  compute  resilient  strains  the  material  properties  as  shown  in  Figure  8  were 
assumed  (from  procedure  given  in  Reference  ?)  and  a  modified  version  of  the 
Chevron  computer  program  was  employed  as  a  mathematical  model  for  computing 
resilient  strains.  From  the  previously  discussed  laboratory  data  the  relation¬ 
ship  of  permanent  to  resilient  strain  for  levels  of  100,  1000,  and  10,000  strain 
repetitions  as  shown  in  Figure  9  was  assumed.  These  relationships  must  be 
considered  as  a  pure  guess  at  the  true  rela'ionships  between  the  computed 
resilient  strains  and  the  resulting  permanent  strains  in  the  pavement  system. 

The  computations  for  determining  the  permanent  strains  are  shown  in  Tables  2, 

3,  and  h,  and  resulting  distributions  of  permanent  strain  with  depth  for 
repetition  levels  of  ]00,  1000,  and  10,000  coverages  are  shown  in  Figures  10, 

II,  and  12  for  pavement  sections  1,  2,  and  3,  respectively.  It  is  to  be 
noted  that  it  has  been  assumed  that,  no  permanent  strain  occurred  in  the 
asphalt  concrete  and  thn*  one  coverage  produces  one  strain  repetition.  It 
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TABLE  B3 
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TABLE  B4 


B  2  2 


PERMANENT  STRAIN  €p 


has  also  been  assumed  that  the  stabilised  gravel  had  the*  same  relationship 
between  permanen.  and  resilient  strain  as  the  crushed  limestone. 

To  compute  the  total  permanent  deformation  the  function  representing 
the  distribution  of  the  permanent  strain  was  numerically  integrated  by  com¬ 
puting  the  area  under  the  distribution  curves.  For  this  computation,  the 
permanent  strain  was  assumed  to  be  zero  at  a  depth  of  2~k0  in.  The  total 
permanent  deformations  for  the  pavement  sections  at  the  different  coverage 
levels  is  given  in  Table  5.  The  comparison  of  the  computed  permanent  deforma¬ 
tion  with  the  measured  permanent  deformation  is  shown  in  Figure  13.  For  each 
of  the  test  sections  the  permanent  deformation  was  overestimated  at  100  cov¬ 
erages  and  underestimated  at  failure. 

There  are  deficiencies  in  the  procedure,  notably  the  assumed  relationships 
between  resilient  strain  and  permanent  strain,  the  inability  of  layered  elastic 
theory  to  predict  strains  and  the  use  of  a  coverage  as  a  strain  repetition, 
which  affected  the  predicted  results.  Of  these  deficiencies  probably  the  most 
serious  and  the  most  difficult  to  correct  is  the  inability  of  the  analytical 
model  to  predict  the  resilient  strains.  The  accuracy  of  the  prediction  is 
greatly  affected  by  the  accuracy  of  the  computed  resilient  strain  which 
emphasizes  the  need  for  more  accurate  models  for  predicting  responses  of 
pavement  systems. 
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TABLE  B5 

Comouted  Surface  Deformation 


Pavement 

Section 

@  100 
Coverages 

@  1000 
Coverages 

§  10,000 
Coverages 

1 

1.1*52 

1.785 

2.31*1* 
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0 . 866 

1.059 

1.285 
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0.772 

0.933 

1.11+3 

SURFACE  DEFORMATION,  IN 


Figure  B13.  Comparison  of  computed  permanent  deformation 
with  measured  permanent  deformation 
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